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The recurrence times between extreme events have been the central point of statistical analyses in 
many different areas of science. Simultaneously, the Poincare recurrence time has been extensively 
used to characterize nonlinear dynamical systems. We compare the main properties of these statis- 
tical methods pointing out their consequences for the recurrence analysis performed in time series. 
In particular, we analyze the dependence of the mean recurrence time and of the recurrence time 
statistics on the probability density function, on the interval whereto the recurrences are observed, 
and on the temporal correlations of time series. In the case of long-term correlations, we verify the 
validity of the stretched exponential distribution, which is uniquely defined by the exponent 7, at the 
same time showing that it is restricted to the class of linear long-term correlated processes. Simple 
transformations are able to modify the correlations of time series leading to stretched exponentials 
recurrence time statistics with different 7, which shows a lack of invariance under the change of 
observables. 
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I. INTRODUCTION 



Recurrence time analysis is a powerful tool to charac- 
terize temporal properties of well defined events 0, • 
It has been recently extensively performed in a rich va- 
riety of experimental time series: records of the cli- 
mate seismic activities @, solar flares 0, spikes 
in neurons , turbulence in magnetic confined plasma 
and stock market indices 01 • Calculated essentially in 
the same way, these analyses receive different names: 
waiting time distribution, interocurrence time statistics, 
distribution of interspike intervals, distribution of lami- 
nar phases, etc. Through a unified perspective, we dis- 
cuss the main properties of these statistical methods, 
which allows us to reinterpret and specify many previ- 
ous results. By a careful discussion of the relevance of 
the probability density function (PDF) of the time series 
data we can easily understand and, in a particular case, 
reject results on earthquake statistics (Sec. [H]). I n the 
case of long-term correlated linear time series we obtain 
a closed expression for the stretched exponential distribu- 
tion of recurrence times which is valid for different recur- 
rence intervals. We show also the lack of invariance of the 
long-term correlations of the time series under transfor- 
mations that simulate the choice of different observables 
of the system (Sec. Illlf) . Before reporting these results 
we start with a proper definition of time series recurrence 
times, we compare it to the Poincare recurrences and we 
discuss briefly some important recent applications of the 
recurrence time analysis. 
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A. Recurrence times 

Time series recurrences 

Assuming a time series point of view, in this paper we 
study the statistical properties of the recurrence time T. 
Given a time series {x n },n = 1, . . . , N, and having de- 
fined a recurrence interval J as a subset of the data range, 
then the ith recurrence time Tj is the time interval An be- 
tween the ith and the i + 1st visit of a time series point 
in /. The recurrence time statistics (RTS) is obtained 
as the distribution P(T) of the sequence of recurrence 
times Tj. 

Evidently, the sequence of recurrence times generated 
this way depends sensitively on the choice of /, which in 
fact will be one prominent issue of this paper. While for 
the recurrence of extreme events the recurrence interval 
is defined by the points above a threshold 

Iext(q) = [Q,°°[, (1) 

in a more general way it may be defined around a posi- 
tion X c with a semi- width S 0, El 

I(X C ,6) = [X C -5,X C + S] . (2) 

Both kinds of intervals are illustrated in Fig. ^ 

Poincare recurrences 

In dynamical systems' theory, another concept of recur- 
rence, known as Poincare recurrences, plays a central 
role. Given a closed Hamiltonian system with ergodic- 
ity on the energy shell, the famous Poincare recurrence 
theorem asserts that almost all trajectories (except for 
a set of zero measure) started inside some subset V of 
the phase space T will return to it infinitely many times. 
In the limit of vanishing volume of this subset, the time 
between consecutive recurrences is the Poincare recur- 
rence time. Despite the well known debates about the 
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FIG. 1: Gaussian distributed time series with recurrence in- 
tervals with measure fj,(I) = 0.01, what implies (T) = 100. 
One extreme interval I ex t(q = 2.3263) — [2.3263, oo[ and 
two recurrence intervals h (-2, 0.0922) = [-1.9078, -2.0922], 
h (0,0.0125) = [-0.0125,0.0125] illustrated in the (a) time 
series and (b) probability density function. 



this paper, is the study of the waiting time between earth- 
quakes or avalanches in models exhibiting self-organized 
criticality (SOC). The idea of studying recurrences in 
SOC started with the first connections between SOC and 
earthquakes 01 • More recently, the investigation of seis- 
mic catalogs of different regions of the globe indicate the 
existence of a universal distribution of recurrence times 
between big earthquakes 6], which may be roughly de- 
scribed as a power-law distribution 



P(T) cx T~ 



(3) 



followed by a faster decay. Simple SOC models have a 
Poisson (exponential) distribution of recurrence times, 
what was used as argument against the use of SOC to 
model not only earthquakes but also (and originally) 
solar flares p|. However, non-Poissonian distributions 
are obtained in more sophisticated SOC models 0,EJ, 
what keeps open the debate over the use of SOC in these 
fields, with the RTS as one of its central ingredients. 



foundations of statistical mechanics (Zermello paradox) 
these ideas motivated throughout the years also mathe- 
matical studies [T^] and, more recently, applications of 
recurrence analysis to many different dynamical systems 
(see Ref. |l3( and references therein). 

Surprisingly, as far as we know, no connection be- 
tween the two recurrence approaches described above 
were made until now. The most evident way to estab- 
lish this relationship is to define an observable x = x(j), 
when rf(t) £ T is the trajectory in phase space of the 
Hamiltonian system. The recurrence volume V is mapped 
to an interval Iy on the real axis by the observation func- 
tion x(ff). However, the sequence of recurrence times of 
the series x n := x("/(t = nAt)) with respect to Iy is 
generally not identical to the sequence of Poincare recur- 
rences of 7(i) with respect of V, since there is usually 
a large set V which also maps to Iy due to the non- 
invertibility of a; (7). Moreover, generally Iy will be of 
the kind of I(X C , 8) rather than I ex t(,q)- 

However, as we will show in this paper, the analogy 
with Poincare recurrences motivates issues related to the 
recurrence times of extreme events which will reveal fun- 
damental insight into their properties. Two main results 
will be the lack of invariance of the RTS under change of 
the observable and the fact that long-term correlations 
are not fully characterized by the autocorrelation func- 
tion. 



B. Earthquakes and SOC models 

The recurrence time between extreme events was re- 
cently used in the analysis of different experimental time 
series 0, 0, IE IE @> @ • One of the most important exam- 
ples of this analysis, which is going to be discussed later in 



C. Long-term correlations and recurrence times 

If time series data {x n } are exponentially (short range) 
correlated, the RTS is well known to be Poissonian, i.e., 
exponential for all T, independent of the choice of I (in 
the limit of small interval — > 0) The same result 
applies to Poincare recurrences (including independence 
of V) if the underlying dynamics is hyperbolic, i.e., in 
well defined mathematical way fully chaotic . Also in 
this case, correlations decay fast. Hence, for systems with 
an exponential decay of correlations, details of defining 
recurrence times and further details of the system are 
irrelevant; instead there exists a unique RTS. 

Many time series data have been found to be long- 
term correlated, i.e., their mean autocorrelation time di- 
verges 0, 0, 0. Typically, this situation is character- 
ized in the time series {x n } (assuming (x n ) — 0) by the 
exponent < 7 C < 1 of the power-law decay of the auto- 
correlation function as a function of the time s 

- N-s 

C x {s) = (xiX i+s ) = — } XiX i+s ~ s" 7c . (4) 

iv — s L — ' 

i=l 

In a recent paper [|, Bunde et al. analyzed the effect of 
long-term correlations on the return periods of extreme 
events, i.e., of recurrence times obtained using recurrence 
intervals of type Q . The main results of Ref. 0, for 
long-term correlated time series can be summarized by 
the following three points. While the first was obtained 
considering statistical arguments, the two others were 
based on numerical simulations. 

(i) The mean recurrence time is equal to the inverse of 
the fraction of extreme points in the series 

N extreme 
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(ii) The statistics of T follows a stretched exponential 

lnP(T) oc -(T/(T)y , (5) 

where 7 = j c is identical to the correlation exponent in 
Eq. |gj). 

(m,) The series of recurrence times is long-term correlated 
with an exponent jt close to j c . 

Statement (i) is the time series analogous of Kac's 
Lemma (Sec.[nJ, statement (ii) will be verified carefully 
(Sec. IIII B|) once we have established the full functional 
form of the stretched exponential (jSJl, and statement (Hi) 
seems not to be generally valid (Sec. 1111 D[) . 

Even if one might argue that based on Ref. [ij the va- 
lidity of Eq. (J3J) is established only for the class of the 
model data chosen there, the reproduction of these find- 
ings for empirical data 0, Q suggests some generality of 
the stretched exponential distribution. Here, the link to 
Poincare recurrences shows the opposite: Hamiltonian 
systems with mixed phase space are long-term correlated 
and show power-law tails in the statistics of Poincare 
recurrence times |13| . In this case, the long-term correla- 
tions are originated by the stickiness of chaotic trajecto- 
ries near the border of integrable islands. They cause a 
kind of intermittent dynamics and manifest themselves in 
complicated higher-order temporal correlations. In fact, 
the temporal properties of typical data are not fully spec- 
ified by the autocorrelation function, Eq. Q, what ex- 
plains why there cannot be a unique RTS for long-term 
correlated data. Connections between the long-term cor- 
relation exponent 7 C and the RTS have to be established 
independently in every class of long-term correlated dy- 
namical systems, as was done for Hamiltonian systems 
with mixed phase space |l3j and fractal renewal point 
processes |2J|. We argue in Sec. IIII 51 that the results of 
Ref. described above are valid for long-term correlated 
linear time series |2| . In this paper we propose a closed 
expression for the RTS of time series of this class, which 
is valid for recurrence intervals of both types and J5J). 

II. MEAN RECURRENCE TIME 

The mean recurrence time 




is a direct result of the choice of the recurrence interval. 
In area preserving dynamical systems Kac's lemma ^1 
states that the inverse of (T) is equal to the ergodic mea- 
sure of the recurrence interval p{I). In the case of sta- 
tionary time series, as illustrated in Fig.^ an equivalent 
result is obtained from the normalized PDF p(x), 

At f 

— = p(I) = / p(x)dx, (6) 
K 1 ) Jxei 

where At is the sampling rate used to record the time 
series [27|. This is the most important constraint to the 



statistics of recurrence times. In the time series analysis 
this measure is estimated as the fraction of valid events 
(points inside the recurrence interval) = N events / N . 
Intuitively, relation © states simply that the total ob- 
servation time t is given by 

t = NAt = N events (T). 

Besides the RTS P(T) the PDF p(x) of the series of 
points itself is typically used to characterize the time se- 
ries. Contrary to other time series analyses (as the de- 
trended fluctuation analysis discussed below), the RTS 
is independent of the PDF. In particular, it is irrelevant 
whether the second moment of the PDF is finite. A time 
series with a well behaved (Gaussian) PDF can have ei- 
ther exponential or power-law RTS 28j. Reversely, a 
time series with fat tails in the PDF can lead to a RTS 
that might be Poisson or power-law [2^| . The reason for 
this is simple: the RTS depends on the sequence of the 
time series points and changes under their temporal rear- 
rangement, which does not change the PDF of the data. 
While the RTS is independent of the PDF of the series, 
the opposite happens to the mean recurrence time (T). 
Once the recurrence interval is defined, whether by re- 
lation (JIJ, J2J or by any other possible definition, the 
PDF p(x) provides (T) through relation JfJJl. 

These two apparently trivial observations, i.e., inde- 
pendence of the RTS and dependence of (T) on the PDF, 
shed new light on previous results. In what follows, we 
exemplify these points in the analysis of recurrence times 
between earthquakes, already mentioned in Scc. lIBI De- 
spite (or because of) the complexity of this field it has 
an important simplicity: the Gutenberg-Richter law 

p(M) oc e- b ln ( 10 > M , (7) 

where A, b are constants and M is the magnitude of the 
earthquake, which is proportional to the logarithm of the 
released energy. The constant b is almost the same for 
different parts of the world and the empirical law J7J is 
valid for 2 < M < 7.5. From our perspective this means 
that the PDF of the time series of seismic activity is 
given [33. 

The mean recurrence time between earthquakes of a 
given magnitude M is obtained inserting the PDF given 
by (7J in relation @), and using the interval of the 
type © with X c = M, 

(T)(M) = To e bln{w)M , (8) 

where To oc (-^-"^lo)^ ■ This relation is equivalent to 
the one obtained previously through a "mean-field ap- 
proach" 0. In Ref. [HI 

it is noted the "remarkable" 
scaling of (T)(M > M c ) oc 10 hMc , which is nothing else 
than a consequence of relation 10 when intervals of the 
type QJ are used with q = M c . 

So far, the relation between (T) and the PDF was used 
to show that the mean waiting time between earthquakes 
is directly related to the Gutenberg-Richter law, but has 
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nothing to do with temporal correlations between earth- 
quakes. On the other hand, the RTS obtained from earth- 
quakes records [|| is an independent result that can be 
used as a test for the dynamical models of earthquakes. 
Recently, it was suggested that in SOC models the se- 
quence of avalanches is uncorrelated 0, E3 (see for 
a counterexample) and should thus be discarded. The 
solution of this debate is beyond the scope of this paper. 
Nevertheless, we note that, as a consequence of the un- 
relatedness of p(x) and P(T), shuffling data of whatever 
distribution randomly (as was done for the time series 
of seismic activity in Ref. 0|) trivially implies P(T) of 
being exponential, also for finite recurrence intervals [llj . 

III. STATISTICS OF RECURRENCE TIMES 




Closed expression of the stretched exponential 
distribution 



We generalize the distribution proposed in Ref. [l| for 
the RTS of long-term correlated time series. Motivated 
by result (it), mentioned in Sec. II Cl suppose that the 
stretched exponential distribution 



P 7 (T) = 



-(6T)i 



(9) 



is valid for all recurrence times T S ]0, oo[. This is 
actually a stronger assumption than Eq. (JSJ. As any 
RTS, Eq. @ must satisfy the following two conditions: 
normalization 

/ P{T)dT = 1, 
Jo 

and the analogous of Kac's lemma © 



poo 1 

(T) = / TP(T)dT = — 

Jo /4 



Imposing these two conditions to the distribution 0, 
it is possible to express the constants a and b as func- 
tions of 7 and fj.(l). Further simplification is obtained 
performing the following transformation of variables r = 
< ^ > = fj,(I)T, i.e., counting the time in units of the 
mean recurrence time. The complete stretched exponen- 
tial distribution for recurrence times is then written as 
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p 1 (r) = a 7 e ^ t)T , with ■ 



' r(i/ 7 )> 
(2 1 /-') 2 r(^i) 

2^ ' 



(10) 



and depends exclusively on the exponent 7. 

Equation (|10|l is illustrated in Fig. [21 for different val- 
ues of 7 in two different ways. Graph (a) (log-log) shows 
that decreasing the value of 7 the distribution starts from 
the exponential (Poisson) case (7 = 1) and approaches 
a power-law (7 — > 0) with an exponent a — 1.5. Graph 



FIG. 2: The stretched exponential distribution 1101 for dif- 
ferent values of the exponent 7. 



(b) shows the distribution in the form that the stretched 
exponentials are seen as straight lines [J Hi- Generally, 
to obtain graph (b) from (a) one needs to divide the dis- 
tribution P(T) by the correct pre-factor a = P(0), which 
is typically unknown. Distribution l|l(J|) shows the de- 
pendence of the pre-factor a on the exponent 7 when the 
stretched exponential p 7 (r) is valid in the whole interval 
of times. For experimental or numerical data, where nei- 
ther a nor 7 are known a priori, the relation between both 
is useful to correctly visualize and fit the RTS. We note 
that in practice the numerical fitting of the exponent 7 
is very sensitive and typically depends on the choice of 
the pre-factor a. 



B. Numerical results for long-term correlated 
linear time series 



We compare now the stretched exponential distribu- 
tion (|10fl to the numerical results of the RTS obtained 
in long-term correlated time series. As in Ref. Q, the 
data were generated using the Fourier transform tech- 
nique 21]: imposing a power-law decay on the Fourier 
spectrum 



f x (k) cx k~P , 



(11) 



with < (3 < 0.5 and choosing phase angles at random, 
we obtain through an inverse Fourier transform the long- 
term correlated time series in x with 7 C = 1 — 2/3 in 
Eq. (@J. The data are Gaussian distributed with (x) = 
0,cr = 1, and Eq. © was used to calculate the times 
r = T/(T). 

Having specified the power spectrum or, correspond- 
ingly, the autocorrelation function for sequences of Gaus- 
sian random numbers means to have fixed all parameters 
of a linear stochastic process. Hence, in principle, the co- 
efficients of an auto regressive (AR(r)) or moving average 
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(MA(r)) process can be uniquely determined, where, due 
to the power-law nature of spectrum and autocorrelation 
function, the orders r of either of these models have to be 
infinite [22]. Hence, the following results are valid for the 
class of linear long-term correlated processes . In other 
words, higher order correlations for this class of processes 
follow trivially from the two-point correlations. 



FIG. 3: (color online) RTS of long-term correlated linear 
time series with N = 2 25 ~ 3 10 7 points for different values 
of 7 C (symbols). Lines are the stretched exponential distri- 
bution 11011 with 7 = 7c- The recurrence interval is extreme 
With fl(Iext) = 10" 2 . 
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We show in Fig. [3] that the stretched exponential dis- 
tribution (|10() with 7 = 7 C describes considerably well 
the RTS, obtained using extreme intervals (Eq. (|TJ), of 
long-term correlated linear time series. The agreement 
is especially good for small values of j c (long correla- 
tions) and q — > oo (which is equivalent to //(/)—> 0). 
This result is a generalization of the result (ii) [ij since, 
using Eq. (|10fl and considering 7 = 7 C , the comparison 
between the theoretical and numerical distributions has 
no free parameter and no fitting is made. 




FIG. 4: (color online) RTS of long-term correlated linear time series with "/ c — 0.1 and different recurrence intervals (centered 
in X c with measure fi(I) = 10" 3 ). The lines are stretched exponentials distributions and the symbols connected by lines are 
the numerical simulations. From (a) to (b) we use the values given by a 7 of the best fitting of Eq. 11011 in (a). In (c) we analyze 
the case X c = for different values of nil), from bottom to top: = lO^ 1 (shifted down by 10 2 ), fi(I) = 10~ 2 (shifted by 

10) and fi(I) = 10 -3 . The gray lines are the Poissonian distribution (7 = 1 in Eq. llOll h a 7 =i = 1 was used for all three cases. 
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Furthermore, we verify in Fig. 0] that, for small 
the distribution is valid also for recurrence intervals 
in the inner part of the data range (centered at X c and 
defined by Eq. @). When X c — > oo, approaching the 
extreme interval, the value of 7 in Eq. (|10|l approaches 
the value of the correlation exponent j c . Decreasing the 
value of X c towards the mean value of the PDF ( (x) = 0) 
results in an increase of 7. This case was analyzed care- 
fully in Fig.^t, where the dependence of the RTS on the 
size of the recurrence interval fi(I) is shown. While for 
big intervals the stretched exponential seems not to hold, 
when — > (the limit Poincare was interested in) the 
distribution for X c — tends to the upper limit 7=1, 
the Poisson distribution. 

In summary, the RTS of long-term correlated linear 
time series with exponent j c , in the limit of small inter- 
val — > 0, is described by the stretched exponential 
distribution (Eq. (|10|) ') for all recurrence times T and for 
recurrence intervals of both types Q and (0). The ex- 
ponent 7 is a continuous and monotonically decreasing 
function of the center X c of the recurrence interval, with 
the limits 



7 



7 C when X c 
1 when X r 



► 00 (extreme), 
0. 



(12) 



This result has a simple interpretation in terms of the 
long-term correlations contained in the time series. Cal- 
culating the RTS to a specific interval measures the cor- 
relation between events inside this interval. In this sense, 
our result suggests that the long-term correlations of the 
time series are concentrated in the extreme events (large 
fluctuations) and vanish for events near the mean value 
(small fluctuations). Relation 11 2[) can then be inter- 
preted as: approaching pure extreme events (/i(/) — * 
and X c — ► 00) the RTS shows the whole correlation and 
thus 7 = 7 C . Approaching pure middle events (/i(J) — > 
and X c = 0) no correlation is detected and consequently 
the Poisson distribution (7 = 1) is recovered. 



C. Change of observables 

The link between recurrence times on time series and 
Poincare recurrences (Sec. II A|) motivates the issue of 
the change of observables. All of the empirical data ex- 
hibiting long-term correlations mentioned before repre- 
sent systems which involve a huge number of degrees of 
freedom. Hence, there is a similarly huge arbitrariness 
in choosing a given observation function ^(7), and the 
natural question is what to expect when we change this 
observation function. 

For instance, the correlations in the weather can be 
studied through records of the daily maximum temper- 
ature or of the daily precipitation 0. For the first ob- 
servable, long-term correlations for times larger than 10 
days were found with an exponent 7 » 0.7 for continental 
stations, independent of the location and of the climatic 



zone of the weather station. On the other hand, the se- 
ries of precipitation, obtained in the same locations and 
for the same time windows, are not long-term correlated. 
A similar situation is observed in financial market data. 
While the fluctuation of prices are typically uncorrelated 
the volatility is long-term correlated 01- This gives al- 
ready a clue that correlations measured on a given time 
series do in fact characterize the fluctuations of the given 
observable but do not characterize the underlying system 
in a more abstract way. 

Here we want to study the dependence of correla- 
tions and RTS on the chosen observable in more de- 
tail by comparing the properties of different observables. 
Generally, both observables x and y are functions of 
the d— dimensional phase space vectors 7, and no sim- 
ple function connecting x and y exists. Since we are 
starting from time series data without underlying multi- 
dimensional phase space, we will restrict the analysis to 
a subclass of changes of observables, where in fact y is 
given by a nonlinear (potentially non-invertible) function 
of x. Hence, we construct time series of different observ- 
ables y as functions of the original long-term correlated 
time series of the variable x. Having in mind a recurrence 
interval defined through (X C ,S) by relation ©, consider 
the following reversible transformation 
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Un 



X n - (X c - 5) 



(13) 



which is essentially the inverse of the original series {x„}. 
If the x-series is Gaussian distributed as considered pre- 
viously, the PDF of the new series {y n } is given by 



p'(y) = 



1 1 



■\{l/y+{Xa-S)) 2 



2tt r 



(14) 



which is illustrated in Fig. for the case X c = 1,5 = 
0.0207. In this figure it is also shown that the interval 1%, 
defined by the same (A c , S) in x, is transformed into an 
extreme interval in y. On the other hand, the extreme 
interval 1% in x is transformed into a recurrence interval 
in the middle of the PDF of y. Since the sequence of 
recurrence times T obtained using the original intervals 
in the x-series is also obtained using the transformed in- 
tervals in the y-series, the RTS remains invariant under 
simultaneous transformation of variables and recurrence 
intervals. Therefore, the previous observation that the 
change of the recurrence interval in the x-series does not 
affect the functional form of the stretched exponential 
distribution (|10fl . but does affect the exponent 7, carries 
over to transformations of the form l|13|) . For instance, 
the RTS of a series obtained from transformation l|13|) 
applied to a time series x with 7 C = 0.1, is well described 
by the stretched exponential distribution (|10f> with (see 
Fig. @J: 7 w 0.55 for extreme interval (I\ in Fig. |SJd) 
and 7 — 0.1 for central interval (I2 in Fig.JSJ). This result 
holds for all reversible transformations. 

An important fundamental question in this context is 
the behavior of the long-term correlations under trans- 
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FIG. 5: PDF of the series (a) of x (Gaussian) and (b) of y 
(Eq. O with X c = 1,8 = 0.0207). The points inside the 
interval Ii(X c = 1,6 = 0.0207) in x become extreme events 
in y. The opposite happens for the extreme interval I2 = 
Jest (« = 2.3263). 



formations of variables. Whereas the normalized auto- 
correlation function remains unchanged under shifts and 
rescalings of x, this is not the case under transformations 
like (|13() . where the transformed time series of y is not 
long-term correlated at all, despite the long-term correla- 
tions of the original x-series (see Fig.[S]where h(2) ~ 0.5). 
We characterize the y-series using the multi-fractal de- 
trended fluctuation analysis [l9j |. which is a much more 
powerful tool than the simple autocorrelation function, 
since for different values of the parameter qDFA different 
scales of fluctuations are amplified. In order to distin- 
guish between the multifractality due to long-term cor- 
relations and due to a broad PDF, the typical procedure 
is to shuffle the time series randomly, i.e., we choose ran- 
domly a new order of the N points of the original time 
series. Since the shuffled series loses all its temporal cor- 
relations but retains the same PDF, the difference be- 
tween the results of the two series (original and shuffled) 
is exclusively due to temporal correlations. In Fig. we 
show the multi- fractal analysis (MF-DFA1 0]) for the 
long-term correlated, Gaussian distributed, linear time 
series {x n } and for the transformed (through Eq. (|13|) ) 
time series {y n }- As expected, in the first case roughly 
a single generalized Hurst exponent h(s) is obtained for 
all scales in the original (h(s) = 1 — 7 c /2 = 0.95) and 
shuffled (h(s) = 0.5) time series. Due to the broad tails 
present in Eq. i|14fl , both the y-series and its shuffled ver- 
sion have multi-fractal spectrum, shown by the nontrivial 
dependence of h(qoFA) on qDFA- The difference between 
the two, which measures the effect of the temporal cor- 
relations, appears for small scales, where the generalized 
Hurst exponent of the shuffled series is smaller. This 
result is consistent with the interpretation made at the 
end of Sec. IIII Bl that the correlations of the x-series is 
concentrated on the extreme events. Through transfor- 
mation (|13|) . the extreme events in x are mapped into 
very small fluctuations in y and the temporal correla- 
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FIG. 6: Generalized Hurst exponent of the time series of x 
and y (N = 2 20 w 10 6 points) as a function of the scale qDFA- 
The horizontal gray lines are the non-correlated value (h = 
0.5) and the expected value for 7c = 0.1 [h = 1 — 7c/2 = 0.95). 
The difference between the original and the shuffled time se- 
ries measures the effect of the correlation at each scale qDFA- 

Through transformation (|13|) we provide an example of 
equivalence between the RTS of different observables ob- 
tained using extreme intervals, and the RTS calculated 
in the same series but using different recurrence inter- 
vals. Always when the transformation of observables is 
invertible, there exist a one-to-one correspondence be- 
tween the original extreme values and a new interval. 
This provides another justification to the generalization 
of the recurrence of extreme events to general recurrence 
intervals, proposed in Sec. II Al inspired by the analogy to 
the Poincare recurrences. 



D. The series of recurrence times 

It is also interesting to apply the distinction between 
the time properties of the series and its PDF, discussed 
in SecTffl to the series of recurrence times {T\, T 2 , ...T/v e } 
itself [25j. In this case this means that the PDF, which 
is the RTS of the original time series, is independent of 
its correlation and shows that the results (ii) and (Hi) 
stated in Sec. II Cl are independent. This is an important 
remark when prediction algorithms are considered, since 
in many cases the correlation between the waiting times 
is more important than their distribution 26] . 

The result (Hi) of Ref. [l| is verified in Fig. [7| through 
the multifractal analysis of the series of recurrence 
times T. Instead of the same correlation exponent we 
find a multifractal spectrum. It is necessarily originated 
by the long-term correlations since the PDF of these se- 
ries are stretched exponential distributions, as verified in 
Fig. which do not have fat tails. 
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FIG. 7: Multifractal spectrum for the series of recurrence 
times ({Ti, Ti, Tjv e }) obtained for intervals with = 
0.01, and different values of X c in long-term correlated linear 
time series with -y c = 0.1 (N e = 2 18 » 2 x 10 5 ). 

IV. DISCUSSION AND CONCLUSION 

Well established regimes of decay of the RTS are ex- 
ponentials and power-laws, and, only recently observed, 
stretched exponentials. We have obtained a closed ex- 
pression for the stretched exponential distribution of re- 
currence times uniquely defined by the exponent 7. As 
limits 7 = 1 and 7 = 0, respectively, we recover the 
exponential and the power-law decay from these (with 
the restriction that the power is fixed to 3/2), suggesting 
that stretched exponentials describe recurrences in sys- 
tems that have neither exponential nor power-law RTS 
but that lie in between these two cases. We have verified 
numerically that the stretched exponential distribution 
is in good agreement with the numerical results obtained 
for long-term correlated linear time series, similarly to 
what was done in Ref. [| . From the point of view of these 
previous results, listed in Sec. II CI we have identified (i) 
with Kac's lemma; generalized (ii) to the stretched expo- 
nential distribution i|l(J|l . which is a function of a single 
parameter and is valid for all recurrence intervals through 
Eq. ljT2|) ; and generalized (Hi), showing that the sequence 
of recurrence times has a multi-fractal spectrum, with an 
exponent 7^ different from j c . In order to verify if the 
fluctuations around the stretched exponential distribu- 
tion, shown in the figures of Sec. lIII 51 are a consequence 
of numerical limitations or real deviations, an analytical 
deduction of the stretched exponential distribution 11 0|l 
is necessary, what remains an open task. 



Performing simple reversible transformations (like 
Eq. 113(1 ) on the original long-term correlated linear time 
series {x n }, we have shown that the stretched exponen- 
tial distribution characterizes also the RTS of extreme 
events in time series that are not long-term correlated. 
The presence and absence of long-term correlations in 
the series of the original observable x and of the trans- 
formed observable y, respectively, is similar to the one re- 
ported above for climatic records (temperature and pre- 
cipitation) and stock-market indexes (volatility and fluc- 
tuation of price). It is remarkable that this interesting 
behavior is obtained already through the simplest possi- 
ble approach, i.e., two different observables that depend 
directly and exclusively on each other. These considera- 
tions emphasize that the temporal characterization of the 
system through the autocorrelation or RTS depend cru- 
cially on the chosen observable. By analyzing both the 
dependence of the exponent 7 of the stretched exponen- 
tial distribution with the center of the recurrence interval 
(relation 1)12(1) and the multi- fractal spectrum of the y- 
series (FIG.|(jJ) we conclude that, in long-term correlated 
linear time series, the correlations are concentrated in the 
extreme events. 

Many interesting questions arise if one supposes that 
the measurements in a given experiment lead to the time 
series of the observable y, introduced in Sec. IIII CI and 
that no natural access to the observable x exists. The 
y-series has a complex multi- fractal spectrum (Fig. EJ, 
a strange PDF (Eq. ((14(1 ) and a non-trivial dependence 
of the RTS with the recurrence interval. Nevertheless, 
through a simple transformation (the inverse of rela- 
tion ((13(1 ) one arrives at the x-series, that has a mono- 
fractal spectrum, is Gaussian distributed and has a sim- 
ple (Eq. ((12(1 ) dependence of the RTS on the recurrence 
interval. This suggests the existence, in some situations, 
of "distinguished observables" where the time series anal- 
ysis is extremely simplified. It is an interesting open 
problem to develop a procedure able to determine the 
transformation (when it exists) that lead to the "distin- 
guished observables" . 
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